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I.  INTRODUCTION 


Recently,  the  question  has  been  raised  as  to  how  the  lightning  threat 
to  aircraft  is  compared  to  that  of  the  high-altitude  nuclear  electromagnetic 
pulse  (NEMP).  The  first  important  step  towards  making  such  a  comparison 
possible  is  to  characterize  the  electromagnetic  properties  of  the  two 
environments.  While  the  high-altitude  NEMP  environment  is  quite  well 
understood,  the  knowledge  of  the  lightning  environment  is  still  in  an 
embryonic  stage.  This  effort  is  to  study  several  possible  models  for 
characterizing  the  lightning  return  stroke  and  to  discuss  how  each  model 
can  be  compared  with  measured  data. 

The  specific  objective  of  this  effort  is  to  calculate  the  natural 
frequencies  of  a  metallic  post  attached  to  a  lightning  channel.  Two 
different  lightning  channel  models  will  be  employed  for  the  calculation. 

The  natural  frequencies  calculated  can  be  compared  with  measured  data 
when  the  latter  become  available.  Such  a  comparison  will  enable  one  to 
determine  the  parameters  inherent  in  the  models.  The  electromagnetic 
characteristics  (such  as  the  fields)  of  a  lightning  channel  can  be  calcu¬ 
lated  and  simulated  with  the  models. 

Section  II  describes  the  two  different  lightning  channel  models. 

Section  III  details  the  approach  which  eventually  leads  to  transcendental 
equations  required  for  numerical  computation.  Explicit  approximate 
analytical  formulas  for  the  natural  frequencies  are  also  derived  in  Section 
III.  The  transcendental  equations  are  described  by  Equations  18  and  19, 
and  the  explicit  analytical  formulas  are  given  by  Equations  20  through  22. 
Section  IV  presents  the  numerical  results.  Section  V  briefly  discusses 
the  application  of  results  to  experimental  measurements. 
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1 1 .  MODELS 


In  this  section  two  lightning  models,  which  may  be  referred  to  as  (1)  the 
resistive  model  and  (2)  the  corona-sheath  model,  are  discussed  for  the 
calculation  of  the  natural  frequencies  of  a  post  with  a  lightning  channel 
attachment.  The  important  features  and  the  required  parameters  of  these 
models  are  described  as  follows. 

1.  RESISTIVE  MODEL 

A  transmission-line  model  has  been  used  extensively  in  the  past  to 
describe  the  electromagnetic  behavior  of  a  lightning  channel  (Ref.  1). 

If  such  a  model  is  made  resistive,  it  would  in  a  sense  account  for  the 
presence  of  a  corona  by  an  introduction  of  a  theoretical  effective  radius. 

The  description  of  a  resistive  model  requires  the  following  two  parameters, 
namely, 

•  the  effective  radius  rg,  and 

•  the  effective  resistance  per  unit  length  R’. 

In  terms  of  these  two  parameters,  the  lightning  channel  together  with  a  corona 
has  a  characteristic  impedance  z£(s)  given  by  (Ref.  2  and  Appendix  A) 

Zj(s)-^*r(s)  (U 


where 


/(s> 


[*"(sre/c)J2 


2  ttR  1  *n(srg/c) 
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s  =  jo)  is  the  complex  frequency,  ZQ,  c,  are,  respectively,  the  free-space 
wave  impedance,  speed  of  light,  permeability,  and  the  superscript  r  is  used 
to  indicate  resistive  model  quantities. 

2.  CORONA-SHEATH  MODEL 

The  corona-sheath  model  has  been  suqaested  in  Reference  3  for  a  lightning 
channel  and  its  surrounding  corona.  In  this  model,  the  lightning  current 


flows  only  in  the  perfectly  conducting  center  channel,  while  all  the  electric 
charges  reside  on  an  effective  corona  surface.  This  model  is  shown  in 
Figure  lb  together  with  a  post.  The  required  parameters  for  the  corona-sheath 
model  are 

•  the  effective  radius  of  the  lightning  center  channel  rQ,  and 

•  the  effective  radius  of  the  corona  surface  rc. 

The  effective  corona  radius  rc  is  related  to  the  charge  per  unit  length  of  the 
corona  Q'  via 


r 


c 


_  IQ'I 

*Vb  * 


Eb  »  3  x  106V/m 


(2) 


where  Eb  and  are  the  air  breakdown  electric  field  and  the  permitti vity  of 
the  surrounding  air.  This  model  can  also  be  considered  as  a  transmission  line 
with  the  characteristic  impedance  (Appendix  A) 

Z‘(s)  *c(s)  (3) 


where 


*c(s) 


and  the  superscript  c  is  used  to  indicate  corona-sheath  model 
quantities.  Clearly,  the  above  two  lightning  models  are  special  cases  of  a 
more  complicated  model  depicted  in  Figure  Ic,  which  allows  for  resistivity  in 
the  center  channel  and,  therefore,  has  the  following  corresponding 
characteristic  impedance  (Appendix  A) 


)- 


2nR'  in  src/c 
*  Tp" 


Z 


/sr^ 


sr 
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Section  III  describes  an  approach  to  solving  this  more  general  model. 
The  two  special  models  given  by  Equations  1  and  3  can  then  be  obtained 
from  Equation  4  by  taking 

•  r0  =  rc  =  re  for  the  resistive  model,  and 

•  R'  *  0  for  the  corona-sheath  model. 
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III.  APPROACH 


The  original  problem  of  Figure  lc  where  the  post  sits  vertically  on  a 
perfectly  conducting  ground  is  equivalent  to  that  of  Figure  2  where  the  images 
are  used  to  account  for  the  ground  effects.  The  natural  modes  of  the  imaged 
problem  can  be  obtained  by  solving  for  the  nontrivial  solutions  of  the 
following  homogeneous  integro-differential  equation  with  certain  appropriate 
boundary  conditions  (Ref.  4  and  the  references  quoted  therein): 


rH  _ 

e-(s/c)Va?  +  (z  -  z')2 


I(z')dz' 


=  0 


z|<  l 


where  a  and  I  are,  respectively,  the  radius  and  length  of  the  post. 
Equation  5  can  first  be  reduced  to  an  integral  equation 


(5) 


e-(s/c)  yja2  +  (z  -  z')2 

. . . .  - - -  -  I(z')dz'  =  B  cosh(sz/c),  (z|  ^  i  (6) 

\l  a2  +  (z  ~z')T 

-l 

where  B  is  the  integration  constant  and  the  following  symmetry  condition 
has  been  imposed 

I(z)  =  I ( ~z )  (7) 

Additional  conditions  are  required  for  solvinq  Equation  6:  they  are  the 
continuation  of  current  and  the  equality  of  potential  at  the  attachment  point 
of  the  post  and  lightning  channel,  i.e.,  at  z  =  l  , 

I_U)  =  I+U)  =  IU) 

♦_(<)  =  ♦+(*)  =  ♦(*)  (8) 

where  the  subscrip  -  or  +  indicates  quantities  at  the  post  side  or  the 
lightning  side. 
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Since  the  potential  at  the  post  side  is  given  by 


*  «  .  ,  1  d_ 

-  4  its  e  dz 

0 


-(s/c)  Va2  +  (z  -  z')2 


Va2  +  (z  -  z'): 


I (z  1  )dz  ' 


(9) 


z  =  i 


and  the  potential  at  the  liqhtning  side  is  equal  to  the  potential  drop  for  a 
current  I+  flowing  through  Zc(s),  i.e.. 


-  I+Zc(s) 


(10) 


Equation  8  is  then  equivalent  to 
rl 


d 
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-(s/c)  Va2  +  (z  -  z*)2 


-l 


Va 2  +  (z  -  z')2 


I(z  1  )dz ' 


4itS  eQZc  (s ) 1  ( 


(ID 


z  =  l 


which  can  be  combined  with  Equation  6  to  solve  for  B  and  give 


el 

-(s/c)  Va2  +  (z  -  z1)2 
5 -  - -  I(z '  )dz' 

Va2  +  (z  -  z')2 

-£ 


sinh(st/c) 


cosh(sz/c ) 


(12) 


This  is  the  final  integral  equation  that  is  needed  to  solve  for  the  natural 
modes  of  the  problem  posed  in  Fiqure  2. 

Since  the  radius  a  of  the  post  is  much  smaller  than  its  length  l ,  an 
asymptotic  theory  can  be  used  to  solve  Equation  12.  First,  rewrite  Equation 
12  in  the  following  form 


12 


1(2) 


dz'  + 


f£I(z •  )e'(s/c)  /2  +  (2  '  z')? 


(z  -  z‘)2 


I(z) 


dz1 


-l 


/a  2  +  (z  -  z(); 


2»(s)I(z) 
sinh ( s  &/ c ) 


cosh($z/c) 


or  equivalently. 


I(z)fi(z) 


+ 


2j£s)lU) 

sinh(st/c) 


cosh(sz/c) 


where 

n(z) 


(13) 


(14) 


(15) 


and  is  a  slowly  varying  function  of  z  with  n(0)  =  2  sinh-1U/a)  -2  tn(2n/a). 
Using  the  fact  that  n(z)  =  n(0)  »  1  and  |i|i(s)|  »  1,  one  can  expand  s  and  I(z) 
in  die  following  fashion 

s  =  s(0)  +s(1)  ♦  . 

I  =  i(0)  +  +  .  (16) 

where  s^m^,  I^171^  are  of  the  order  of  l/d71  or  l/^m.  Substituting  Equation  16 
in  Equation  14,  one  has 


I'°>(z)  =  A 

(0) 


Q(0) 
o  b(z) 


[-  cosh  z/c j 


1  +  — -  cothfs^  i/c]  =  0 


and 


(17) 

(18) 
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nfol  +  c  co'h^ 

d’i-U1 

0 

>) 

i(U)  w(  l)  l  cosh 

V  c  ) 

B 

SSH 

— 

where  (Ref.  5) 


E(  c)  =  Ejfc)  +  in  ?  +  t 


^U)  =  j  dt  =  exponential  integral 


Y  =  0.577....  =  Euler's  constant 


When  |iKs)|  Equations  16  throuqh  19  aqree,  as  they  should,  with 
Equations  16  and  17  of  Reference  4  where  a  post  without  a  liqhtninq  attachment 
was  discussed.  Equation  18  can  be  used  to  solve  for  s^,  from  which  s^  can 
then  be  calculated  via  Eauation  19. 

When  R'  i  »  ZQ  and/or  a  »  (rQ,  rc ,  rg),  one  has  |i/>(s)|  »  fi(Jt),  and 
Equation  18  can  be  approximately  solved.  The  approximate  solution  is 


s(0)* 

J] _ j(2n  +  1 )  -m 

c  2 


>  n  =  0,1,2, 


where  yp  =  <Ks)  with  s  =  ;'(2n  +  1 )  ttc / ( 2 z) ,  i.e.. 


*n  21  j  (*»  [t?TT1T^]  ■ j?)  (*  [(Tn'4r)?n;]  -*?) 

R'4  4.1  Lj  2 1  1  ■ 

'  T~  (2n"+n  r1 


Furthermore,  Equation  19  can  be  approximated  by 


1 

=  •  TaflJT 


JZn  [( 2n  +  l)2ir] 


+ 


(22) 


In  Equations  20  through  22,ft(t)  =  tn(4t/a),  0.(0)  ~2£n(24/a),  and  the 
subscript  n  indicates  the  quantity  referring  to  the  n-th  natural  mode. 

The  sum,  s^  +  s^,  obtained  from  Equations  20  through  22  gives  the 

approximate  analytical  formula  for  the  natural  frequencies  of  a  post  with  c 

lightning  channel  attachment.  The  accuracy  of  the  approximate  analytical 

formula  will  be  checked  against  the  numerical  results  computed  from  Equations 

18  and  19  in  Section  IV.  It  is  seen  from  Equation  20  and  22  that  the  effect 

of  a  channel  attachment  on  the  post's  natural  frequencies  is  entirely  accounted 

for  by  the  last  term  of  Equation  20.  This  term  is  proportional 
-1/2 

to  (R'Ji/Z  )  when  the  lightninq  channel  is  highly  resistive,  and  vanishes, 
as  expected,  when  R'  *  <*>. 


IV.  NUMERICAL  RESULTS 


m 


Equations  1A  and  19  were  numerically  solved  on  a  microcomputer  for 
4<( s)  qiven  by  Equation  4.  The  results  of  s^  +  s^  are  given  in  Figures  3 
through  14  for  post  radius  a  =  0.05  m,  and  various  post  length  l  and  lightning 
parameter  values.  The  lightning  parameters  are  the  corona  radius  rc,  the 
channel  radius  rQ,  the  combined  effective  radius  rg,  and  the  effective 
resistance  R'  per  unit  length  (Fig.  1). 

Figures  3  through  9  are  the  natural -frequency  plots  for  the  resistive 
model  with  r0  s  rc  =  re  and  nonzero  R'£/Zo.  Figures  3  and  4  show  the  effect 
of  the  channel  resistance  R'  on  the  post's  natural  frequencies  for  two 
different  values  of  effective  lightning  radius  rg.  The  imaginary  part  of  the 
natural  frequencies  remains  almost  unchanged  as  R'z/Zo  varies.  On  the  other 
hand,  the  absolute  value  of  the  real  part  (i.e.,  the  decaying  constant) 
increases  when  R'X./Zo  decreases.  This  is  expected  because  a  post  with  a 
lightning  attachment,  of  lower  resistance  is  different  from  a  post  without  a 
lightning  attachment. 

Figures  5  and  5  show  the  effect  of  the  effective  lightning  radius  ra  on 

a  e 

the  post's  natural  frequencies  for  R '  VZQ  >  10  and  equal  to  10.  The  natural 
frequencies  are  found  to  depend  only  weakly  on  rg.  This  is  because  the 
dependence  of  the  characteristic  impedance  Zc  of  a  lightning  channel  on  rg  is 
logarithmic.  Also,  when  the  channel  is  very  resistive,  i.e.,  when 

4 

R't/ZQ  >  10  ,  the  natural  frequencies  become  independent  of  rg,  and  are  the 
same  as  those  of  a  post  without  a  lightning  attachment. 

Figures  7  through  9  show  the  effect  of  the  post  length  t  on  the  natural 
frequencies.  When  R'  i'lQ  is  large  (i.e.,  >  10^),  the  decaying  constant 
increases  as  i  decreases.  The  reason  is  that  there  is  more  radiation  loss  for 
a  thicker  post  of  a  given  length.  When  R'z/Zo  is  small  (<  10),  the  dependence 
of  the  normalized  natural  frequency  on  l  is  weak.  This  is  attributable  to  the 
fact  that  the  damping  of  the  post  current  is  mainly  caused  by  the  leakage  into 
the  lightning  channel  whose  dependence  on  l  is  logarithmic. 

Figures  10  through  14  give  the  natural  frequencies  of  the  corona-sheath 
model  with  R't/ZQ  =  0.  Figures  10  and  11  show  the  effect  of  the  corona  radius 
rc  whereas  Figures  12  through  14  show  the  effect  of  the  post  length  i  on  the 
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RESISTIVE  MODEL  CALCULATION 


ed  natural  frequencies  of  a  post  of  10m  long,  5cm  in  radius  with 
ina  attachment  of  r  =1  irniand  various  R' .  The  branch  of  poles  with 
10^  also  correspond  to  the  poles  of  a  post  in  free  space. 


numerically  calculated 
results 


Normalized  natural  frequencies  of  a  post  of  10  m  long,  5  cm  in  radius,  with 
a  lightning  attachment  of  R ’ A/Z  =  10  and  various  r  .  In  the  figure,  the 
natural  frequencies  for  a  post  with  no  lightning  attachment  are  also  included 


RESISTIVE  MODEL  CALCULATION 


Figure  7.  Normalized  natural  frequencies  of  a  post  of  various  length  and  5  cm  In 

radius,  with  a  lightning  attachment  of  r  =1  mm  and  R’t/Z  >  10-4.  These  poles 
also  correspond  to  those  of  a  post  in  frle  space.  0 


numerically  calculated  results 


ed  natural  frequencies  of  a  post  of  various  length  and  5  cm  In  radius 
ightning  attachment  of  re  =  1  mm  and  R'£/Z0  =  10.  In  the  figure,  the 
frequencies  for  posts  with  no  lightning  attachment  are  also  included. 
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Figure  9.  Normalized  natural  frequencies  of  a  post  of  various  length  and  5  cm  In  radius 
with  a  lightning  attachment  of  r  =  1  cm  and  R'£/Z0  =  10.  In  the  figure,  the 
natural  frequencies  for  posts  with  no  lightning  attachment  are  also  included. 
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Normalized  natural  frequencies  of  a  Dost  of  10m  long  and  5cm  in  radius  witn  a 
lightning  attachment  of  R'H/2  ~  10"S  r  =  lrrm  and  various  r  .  In  the  figure, 
the  natural  frequencies  for  a  post  with  no  lightning  attachment  are  also  included 


lalized  natural  frequencies  of  a  post  of  10  m  long  and  5  cm  In  radius,  with  a 
tning  attachment  of  small  R'A/Z  ,  r  =  1  cm  and  various  rc.  In  the  figure, 
natural  frequencies  for  a  post  with  no  lightning  attachment  are  also  included. 
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Normalized  natural  frequencies  of  a  post  of  various  length  and  5  cm  in  radius,  with 
a  lightning  attachment  of  small  R'SL/Z  ,  r  =  1  mm  ,  and  r<-  =  1  cm.  In  the  figure, 
the  natural  frequencies  for  posts  witn  no  lightning  attachment  are  also  included. 
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Figure  14.  Normalized  natural  frequencies  of  a  post  of  various  length  and  5  cm  in  radius,  with 
a  lightning  attachment  of  small  R'l/l  and  r  =  r  =  1  cm.  In  the  figure,  the 
natural  frequencies  for  posts  with  no  lightning  attachment  are  also  included. 


post's  natural  frequencies.  The  normalized  natural  frequencies  are  nearly 
independent  of  l. 

The  approximate  analytical  results  given  by  Equations  20  and  22  are  also 
plotted  in  Figures  3  through  14  for  comparison  with  the  numerical  results. 

The  agreement  is  generally  excellent  for  the  lower-order  modes. 
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V.  SUMMARY 


In  /,  l.hr,  r  t.  gi  t.h«-  •  Imj  I  -it  ion  of  the  natural  frequencies 

of  a  metallic  post  on  a  qround  plane  with  a  hohtninq  chanriel  attached  and 
with  various  post  radii  and  post  lenqths. 

As  measured  frequency  data  become  available,  they  can  be  compared  to  the 
model  calculation  presented  in  this  report.  Such  a  comparison  will  enable  one 
to  determine  the  model  parameters,  such  as  the  effective  resistance  per  unit 
lenqth  R‘,  the  effective  radius  of  the  liqhtmng  center  channel  r  and  the 
effective  radius  of  the  corona  surface  rc. 

A  closer  observation  of  all  the  fiqures  presented  in  Section  IV  indicates 
that  there  is  a  definite  relationship  between  the  natural  frequencies  (snt/c) 
and  R',  rQ,  rc.  In  particular,  the  natural  frequencies  move  away  from  the 
jto-axis  when  R'  becomes  greater.  This  trend  enables  one  to  quickly  determine 
the  rouqh  ranqes  of  the  model  parameters  R',  r  ,  rc  for  a  lightmnq  attached 
to  a  post  when  one  comoares  the  poles  extracted  from  the  measured  post  current 
with  these  calculated  poles.  For  example,  if  the  post  lenqth  i  -  10  m,  and 
the  post  radius  a  *  5  cm,  and  if  the  measured  data  satisfy  the  inequality 
I Re( sn£/c ) )  <  0.4,  one  could  quickly  obtain  the  conditions  from  the  trend 
that  r  *  r  <  0.2  cm  and  R '  t/7  >  100,  whereas  the  more  exact  values  for  r  , 

v  U  U 

rc,  and  R'  have  to  be  determined  by  a  comparison  of  measured  poles  with  the 

calculated  poles  such  as  those  in  Fiqure  3.  On  the  other  hand,  if  the 

measured  data  show  |Re(s  */c)|  >  0.4,  the  condition  that  R't/Z  <  1  would  be 

n  0  ~ 

obtained,  while  the  more  exact  values  for  rQ,  rc  and  R'  have  to  be  determined 
by  comparing  the  measured  poles  with  the  calculated  curves  such  as  those  of 
Fiqures  3,  4,  10  and  11. 

From  the  obtained  model  parameters  R',  rQ,  rc  and  the  measured  post 
current,  one  can  synthesize  the  lightning  current  distribution,  from  which 
other  electromagnetic  characteristics  such  as  the  fields  generated  can  be 
calculated. 

I 

The  model  employed  in  this  report  involves  the  assumption  that  the  qround 
is  perfectly  conducting  and  infinite  in  extent.  This  assumption  is  not 
expected  to  introduce  significant  error,  since  the  fields  are  quite  localized 
along  the  lightning  channel  and  the  measuring  post. 
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APPENDIX  A 


CHARACTERISTIC  IMPEDANCE  OF  A  LIGHTNING  MODEL 

In  this  appendix,  the  characteristic  impedance  for  the  lightning  model 
depicted  in  Figure  lc  will  be  derived.  This  impedance  is  reducible  to  the  two 
simpler  cases  corresponding  to  the  lightning  models  of  Figures  la  and  lb. 

In  Reference  2,  it  is  shown  that  a  perfectly  conducting  wire  of  radius  r 
has  a  characteristic  impedance  2°  given  by 

*^Ko(sr/c)  (A.l) 


Since  the  inductance  per  unit  length  L*  and  the  capacitance  per  unit  length  C' 
of  the  wire  is  related  to  Z°  via 


Z°  =  /UTV  =  cL'  *  (cCT1 

v 


(A. 2) 


one  has 


l'  >^K0(sr/c) 

C'  =»  2TreQ/K0(sr/c) 


(A. 3) 


With  tht  so  equation';,  the  charac  ter  i  s  t  ic  impedance  tor  the  lightning 
model  of  Figure  lc.  can  be  formulated.  Since  for  this  model  the  currents 
flow  mainly  in  the  lightning  channel  of  radius  rQ  and  the  charges  reside 
only  on  the  effective  corona  surface  of  radius  rc>  this  lightning  model 
will  have  a  series  impedance  per  unit  length  V  and  a  shunt  admittance  per 
unit  length  Y’  given  by 


